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Abstract 

Hierarchical  nanostructures,  ranging  through  atomistic,  molecular  and  macroscopic  scales, 
represent  universal  features  of  biological  protein  materials.  Here  we  show  for  the  case  of 
alpha-helical  (AH)  protein  domains  that  this  use  of  molecular  hierarchies  within  the  structural 
arrangement  leads  to  an  extended  physical  dimension  in  the  material  design  space  that  resolves 
the  conflict  between  disparate  material  properties  such  as  strength  and  robustness,  a  limitation 
faced  by  many  synthetic  materials.  An  optimal  combination  of  redundancies  at  different 
hierarchical  levels  enables  superior  mechanical  performance  without  additional  material  use. 
Our  analysis  is  facilitated  by  the  application  of  a  Hierarchical  Bell  model  (HBM),  which 
explicitly  considers  the  hierarchical  architecture  of  H-bonds  within  the  protein  structure, 
providing  a  structure-property  relationship  of  strength  properties  of  AH  protein  nanostructures. 
The  HBM  is  validated  by  large-scale  molecular  dynamics  simulations  of  several  model  protein 
structures.  Our  findings  may  enable  the  development  of  self-assembled  de  novo  bioinspired 
nanomaterials  based  on  peptide  and  protein  building  blocks,  and  could  help  in  elucidating  the 
mechanistic  role  of  AHs  in  cell  signaling  and  mechanotransduction. 

Supplementary  data  are  available  from  stacks.iop.org/Nano/20/075103 
(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

The  origin  of  how  naturally  occurring  biological  protein 
materials  (e.g.  spider  silk,  bone,  tendon,  skin)  are  capable 
of  unifying  disparate  mechanical  properties  such  as  strength 
(ability  to  sustain  large  stresses  without  fracture)  and 
robustness  (ability  to  undergo  deformation  without  fracture, 
despite  the  presence  of  defects,  equivalent  to  toughness), 
as  well  as  other  physical  properties  such  as  self-healing 
ability,  adaptability,  changeability,  and  evolvability  into  multi¬ 
functional  materials  is  of  significant  interest.  However, 
the  molecular  basis  of  these  properties  remains  largely 
unknown  [1-3].  Many  synthetic  materials  are  not  capable 
of  unifying  strength  and  robustness,  being  either  extremely 


strong  with  little  ductility  (e.g.  ceramics,  glass,  silicon),  or 
weak  with  extreme  ductility  (e.g.  soft  metals  like  copper)  [4]. 
The  combination  of  these  disparate  properties  into  synthetic 
materials  remains  an  open  challenge  on  the  way  towards  the 
development  of  biomimetic  structures  and  material  designed 
from  the  nanoscale  up. 

The  folded  structure  of  proteins  is  stabilized  by  a  variety  of 
chemical  driving  forces  including  hydrophobic  effects,  H-bond 
formation  as  well  as  charge  interactions  [5].  Simultaneously, 
in  a  folded  state  a  loss  in  configurational  entropy  appears  due 
to  the  loss  of  degrees  of  freedom,  as  almost  all  rotomeric 
angles  of  the  backbone  and  the  side  chains  are  restricted  to  one 
position.  Thus  the  enthalpic  effects  of  the  folding  process  need 
to  be  higher  than  the  entropic  contributions  at  given  conditions 
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(e.g.  temperature,  pH,  etc).  In  alpha-helical  (AH)  based  protein 
structures  (including  individual  AHs,  alpha-helical  coiled-coils 
(CCs),  higher  order  filament  assemblies,  etc)  H-bonds,  ionic 
as  well  as  hydrophobic  interactions  play  an  important  role  in 
creating  and  stabilizing  the  protein  structure.  Whereas  AHs  are 
stabilized  primarily  by  H-bonds,  the  structure  of  CCs  is  related 
to  additional  ionic  and  hydrophobic  interactions  [6,  7]. 

An  AH  is  generated  when  a  single  polypeptide  chain 
twists  around  on  itself,  stabilized  by  H-bonds  made  between 
every  fourth  residue,  linking  the  O  backbone  atom  (hydrogen 
acceptor,  providing  free  electrons)  of  peptide  i  to  the  N 
backbone  atom  (hydrogen  donor)  of  peptide  /  +  4  in  the 
polypeptide  chain.  Consequently,  at  each  convolution,  3.5  H- 
bonds  are  found  in  a  parallel  arrangement  that  stabilize  the 
helical  configuration  of  AHs  [8].  Due  to  the  stabilizing  role 
of  H-bonds  [5],  rupture  of  these  protein  structures  itself  is 
mainly  determined  by  breaking  of  these  H-bonds.  Once  these 
H-bonds  are  broken,  the  protein  has  the  degrees  of  freedom 
necessary  for  further  protein  unfolding.  This  is  in  particular 
the  case  when  the  breaking  of  H-bonds  appears  in  a  highly 
driven  system,  far  away  from  equilibrium.  In  contrast,  at 
very  small  deformation  rates  near  the  equilibrium  the  H-bond 
rupture  needs  significantly  less  force  than  the  protein  unfolding 
itself,  which  is  determined  by  entropic  effects  of  stretching 
the  protein  backbone  (and  restricting  the  degrees  of  freedom 
again)  (see,  e.g.  the  discussion  recently  reported  in  [9]).  In 
this  paper  we  focus  on  protein  rupture  mechanics  in  a  highly 
driven  system,  and  therefore  consider  H-bonded  interactions 
as  the  main  contributor  to  determine  their  strength  properties. 

Although  it  is  recognized  that  H-bonds  play  a  crucial 
role  in  defining  the  strength  properties  of  fundamental  protein 
constituents  in  biological  protein  materials  [2,  7,  8,  10], 
the  role  of  their  characteristic  nano  structured  hierarchical 
arrangement  as  well  as  their  influence  on  the  larger- scale 
mechanical  properties  of  protein  filaments  (that  is,  elasticity, 
fracture/rupture,  energy  dissipation,  etc)  remains  largely 
unknown.  However,  the  molecular  and  supermolecular, 
mesoscale  behavior  of  basic  protein  constituents  is  elementary 
in  order  to  progress  towards  an  integrated  understanding  of 
the  mechanical  role  proteins  play  in  biological  systems,  for 
chemomechanical  coupling,  signaling  cascades,  and  protein- 
protein  interactions  [3,  10-12],  as  well  as  for  the  use  of  protein 
building  blocks  and  hierarchical  structures  in  the  design  of  new 
nanomaterials  [13]. 

A  promising  path  to  build  de  novo  macroscopic  structures 
for  technological  applications  from  nanoscopic  elementary 
buildings  blocks  is  to  arrange  them  systematically  in  a 
hierarchical  order  as  observed  in  Nature.  Thereby  protein 
domains,  macromolecules,  nano  wires  or  different  types  of 
nanotubes  could  serve  as  possible  elementary,  universal 
building  blocks.  However,  engineering  these  structures  and 
controlling  their  chemomechanical  properties  across  multiple 
length- scales  requires  a  fundamental  understanding  of  the 
basic  materials  science  concepts  found  in  these  materials. 
Theoretical  efforts  to  elucidate  these  concepts  are  reported  in 
this  paper. 

Here  we  focus  on  an  analysis  of  simple  alpha-helical 
(AH)  protein  domains  as  model  systems,  which  are  universally 


found  nanostructural  components  of  many  biological  protein 
materials.  These  protein  domains  play  a  crucial  role  in  the 
signaling  and  deformation  behavior  of  cytoskeletal  protein 
networks  in  cells  (e.g.  intermediate  filaments  vimentin  and 
lamin  as  well  as  actin  [7,  8,  11]),  and  in  determining  the 
mechanical  properties  of  hair,  hoof,  feather  and  many  other 
important  structural  protein  materials  [8].  Nanostructured 
AH-based  protein  domains  universally  define  the  nanoscale 
architecture  of  these  protein  materials.  Through  the  analysis 
reported  here,  we  illustrate  that  the  nanostructural  arrangement 
of  universal  protein  building  blocks  plays  a  crucial  role  in 
defining  their  material  performance.  Thereby,  the  occurrence 
of  material  hierarchies  (that  is,  the  arrangement  of  subunits 
to  units,  which  themselves  form  larger- scale  structures,  etc) 
is  a  particularly  important  aspect.  The  analysis  reported 
here  is  focused  on  AH-based  protein  filaments  as  a  model 
system  to  illustrate  fundamental  material  concepts  that  could 
be  translated  to  the  design  of  other  nanostructures  such  as 
hierarchically  arranged  carbon  nanotubes  or  nano  wire  bundles. 

1.1.  Hierarchical  H-bond  structures 

In  biological  tissues,  macroscopically  applied  stress  is 
forwarded  to  microscopic  hierarchical  scales,  where  rupture 
of  H-bonds  mediates  deformation,  thereby  controlling  the 
response  at  the  macroscopic  protein  filament  level.  Since 
rupture  of  individual  H-bonds  equals  a  chemical  reaction 
(breaking  of  the  backbone  H-bond  and  building  of  H- 
bonds  with  surrounding  water  molecules  governed  by  a 
difference  in  H-bond  building  energy,  see  figure  1(c)),  an 
integrated  chemomechanical  approach  is  compulsory  for 
the  understanding  of  protein  fracture  mechanisms  and  the 
development  of  constitutive  mathematical  relations. 

The  lowest  hierarchy  of  AHs  (and  other  protein  structures) 
is  typically  composed  of  arrangements  of  weak  H-bond 
interactions,  organizing  amino  acids  in  stable  elementary 
building  blocks  such  as  AHs,  which  form  hierarchical 
arrangements  such  as  CCs,  supercoils  and  filamentous 
structures  [7,  8,  11,  12]  (see  figures  1(a)  and  (b)  and  2(a), 
exemplified  for  AH-based  structures  as  they  appear  in 
intermediate  filaments  [7]).  The  key  to  enable  the  development 
of  an  accurate  bottom-up  mechanistic  understanding  of  the 
strength  properties  of  such  AH  protein  domains  and  assemblies 
thereof  must  therefore  include  an  explicit  description  of 
nanopatterned  H-bond  arrangements,  by  adapting  a  system 
view  of  materials,  in  the  spirit  of  a  merger  of  nanostructure 
and  materials.  The  concept  of  hierarchical  arrangements, 
where  each  element  consists  of  several  subelements  provides  a 
broadly  applicable,  yet  simple  model  to  describe  the  geometry 
of  AH  protein  structures  [14].  The  schematic  representation 
depicted  in  figures  1(a)  and  (b)  and  2(a)  could  in  principle  be 
used  to  represent  other  structural  proteins  such  as  beta- sheets, 
beta-helices  or  tropocollagen,  since  the  difference  between 
many  protein  structures  is  the  geometrical  arrangement  of  H- 
bonds. 

H-bond  rupture  mechanisms  can  be  seen  as  analogs  to 
the  nucleation  of  dislocations  in  ductile  materials  or  the 
rupture  of  covalent  bonds  in  brittle  materials,  representing 
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Figure  1.  Illustration  of  structural  hierarchies  and  their  representation  in  the  Hierarchical  Bell  model  (subplot  (a)),  as  well  as  representation  of 
the  corresponding  physical  system  (subplot  (b)).  The  inlay  in  the  upper  part  of  subplot  (a)  shows  a  single  AH  structure  with  ^3  H-bonds  per 
convolution.  The  Hierarchical  Bell  model  enables  one  to  predict  the  strength  of  different  hierarchical  bond  arrangements  as  a  function  of  the 
deformation  speed.  Subplot  (b)  shows  the  physical  system  that  is  represented  in  the  hierarchical  model  in  subplot  (a).  Subplot  (c):  statistical 
theory  to  predict  the  bond  rupture  mechanics  via  the  bond  energy  landscape  [15,  17].  The  graph  depicts  the  energy  as  a  function  of 
deformation  along  a  deformation  variable,  along  a  particular  pathway  that  leads  to  bond  rupture.  Here  /  is  the  applied  force,  and  Xb  is  the 
displacement  in  the  direction  of  the  applied  force,  corresponding  to  the  lateral  displacement  that  is  necessary  to  overcome  the  bond  breaking 
distance  of  a  H-bond.  Given  that  is  the  distance  to  break  a  single  H-bond,  the  distance  =  XbCos  0  denotes  the  lateral  displacement  at 
bond  breaking,  with  the  angle  0  as  the  angle  between  pulling  direction  and  orientation  of  the  H-bond  inside  the  molecule.  This  fundamental 
view  of  single  bond  behavior  is  scaled  up  through  several  hierarchies  in  the  Hierarchical  Bell  model. 


fundamental  unit  deformation  events  [4,  15-19].  Notably, 
larger- scale,  effective  properties  of  these  hierarchical  H-bond 
structures  can  not  be  calculated  by  conventional  mean-field 
averaging  approaches,  not  only  due  to  an  insufficient  number 
of  subelements  [14],  but  also  since  information  may  be 
forfeited  that  is  crucial  for  the  structure’s  behavior  several 
scales  up.  Recent  results  provide  strong  evidence  that  the  key 
to  understand  the  mechanical  response  of  AH-based  protein 
structures  is  to  consider  the  rupture  dynamics  of  H-bonds  at 
mesoscale  [17,  18,  20-25]. 

1.2.  Outline  of  this  paper 

The  central  question  addressed  in  this  paper  is,  how  does  a 
hierarchical  AH  protein  structure  respond  to  mechanical  load 
and  how  does  this  response  relate  to  the  protein’s  structure? 
Further,  how  does  the  hierarchical  arrangement  control  the 
multi-scale  process  of  fracture  and  deformation  in  protein 
materials?  How  can  AH  protein  domains  overcome  the 
competition  between  strength  and  robustness?  In  this  paper  we 
address  this  question  by  developing  a  simple  theoretical  model 
coupled  with  molecular  dynamics  (MD)  simulations.  We  begin 
with  the  presentation  of  the  theoretical  model  and  validation 
simulations.  We  then  continue  with  a  discussion  of  several 
case  studies,  all  focused  on  illustrating  the  use  of  hierarchies 
in  achieving  improved  material  properties  that  exceed  those  of 
the  elements  alone.  We  conclude  the  paper  with  a  discussion 
of  implications  for  bioinspired  approaches  in  nanotechnology. 


2.  Theoretical  and  computational  modeling 

In  this  section  we  summarize  theoretical  and  computational 
methods  used  for  the  analysis  reported  in  this  paper. 

2.1.  Hierarchical  Bell  model 

The  Bell  theory  [26]  has  been  used  successfully  to  describe 
the  rupture  mechanics  of  adhesion  bonds,  but  has  also  been 
used  to  describe  the  rupture  dynamics  of  H-bonds.  The  central 
element  of  Bell’s  theory  is  the  off  rate,  which  describes  how 
often  a  bond  dissociates  per  unit  time  (this  concept  is  illustrated 
in  figure  1(c)): 

/  (fib  -  /XbCOS(0))\ 

- V — ) 

In  order  to  allow  capturing  the  pulling  speed  dependence 
of  the  strength  of  H-bonds,  here  we  link  the  off  rate 
to  the  bond  breaking  speed,  by  multiplying  the  off  rate 
(reciprocal  of  the  bond  breaking  time)  by  the  distance  Xb, 
which  needs  to  be  overcome  in  order  to  break  the  bond: 
V  =  x-^b-  Despite  the  usefulness  of  this  approach, 

existing  models  derived  on  Bell’s  concept  are  only  capable 
of  treating  conglomerates  of  bonds,  as  they  only  consider 
an  effective  energy  barrier  that  corresponds  to  multiple 
individual  bonds.  Moreover,  thus  far  all  attempts  have  failed 
to  predict  the  strength  of  hierarchical  arrangements  of  H- 
bonds  from  fundamental  principles,  explicitly  considering  the 
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Figure  2.  Validation  of  Hierarchical  Bell  model  via  direct  atomistic 
simulation  in  explicit  solvent.  Subplot  (a):  overview  of  three  model 
systems,  including  single  AH,  coiled-coil  with  two  AHs  (CC2),  and  a 
coiled-coil  protein  with  four  AHs  (CC4).  The  schematic  on  the  right 
shows  the  hierarchical  representation  of  each  physical  system 
(following  the  approach  shown  in  figure  1).  Each  blue  line  represents 
a  single  H-bond  (in  the  atomistic  simulations,  all  structures  are 
embedded  entirely  in  a  water  skin  during  deformation;  however, 
water  molecules  are  not  shown  here  for  clarity).  Subplots  (b)  and  (c): 
force-extension  curves  and  unfolding  force  as  a  function  of  pulling 
speed.  The  beginning  of  the  plateau  regime  (regime  II),  following  the 
initial  elastic  stretching  (regime  I)  in  subplot  (b)  defines  the  unfolding 
force,  the  quantity  plotted  in  subplot  (c).  The  continuous  lines  in 
subplots  (c)  are  the  predictions  from  the  Hierarchical  Bell  model, 
equation  (2).  All  predictions  have  been  made  based  on  the  same 
input  parameters  of  H-bond  energy  and  therefore,  the  differences  are 
solely  due  to  hierarchical  effects.  It  is  noted  that  the  force  level  in  the 
plateau  regime  is  defined  by  the  unfolding  force,  and  therefore,  the 
dissipated  energy  is  proportional  to  this  force  value. 


geometrical  arrangements  of  the  bonds.  Here,  we  introduce 
the  development  of  the  Hierarchical  Bell  model  that  overcomes 
these  limitations.  A  detailed  derivation  of  the  mathematical 


model  is  included  in  the  supplementary  material  (available  at 
stacks.iop.org/Nano/20/075 103). 

As  discussed  above,  the  ‘original’  Bell  model  (equa¬ 
tion  (1))  does  not  distinguish  between  a  single  chemical  bond 
and  protein  architectures  that  include  several  fundamental 
bonds.  For  instance,  whether  a  single  H-bond  ruptures  or  if 
several  H-bonds  rupture  simultaneously  is  captured  in  an  ef¬ 
fective  value  of  E\y  \  however,  this  change  in  mechanism  is  not 
explicitly  noted  in  the  model  given  in  equation  (1)  and  subse¬ 
quent  expressions.  In  order  to  estimate  the  strength  and  the 
energy  landscape  of  a  protein  without  performing  any  simu¬ 
lations  or  experiments  and  thus  to  make  the  model  predictive, 
the  model  is  extended  to  explicitly  consider  the  structural  hier¬ 
archies  of  the  protein  structure  with  the  only  input  parameters 
being  the  energy  of  a  H-bond  and  the  rupture  distance.  The 
H-bond  represents  a  hierarchical  structure,  reaching  from  indi¬ 
vidual  H-bonds  at  the  lowest,  atomistic  level  to  a  collection  of 
H-bonds  at  the  next  higher,  molecular  protein  scale. 

In  this  model,  the  system  breaking  force  fhn  and  the 
energy  barrier  of  a  system  consisting  of  n  hierarchies  (the 
mathematical  symbols  are  explained  in  table  1)  is  given  by 


fhn  = 


k^T 

Ab  COS  6 


X  In 


bi-i 

ki-i 
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Ab  COS  6 
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[5(‘; 


In 


bi-x 

ki-i 


+  ln 
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These  equations  now  enable  us  to  predict  the  unfolding 
force  at  any  pulling  speed,  once  the  structural  geometry  and 
the  energy  landscape  of  a  single  H-bond  is  known. 

We  note  that  this  model  only  considers  H-bonds  as 
structural  elements  in  the  definition  of  AH-based  protein 
structures,  representing  a  limitation  of  this  model.  Thereby 
it  does  not  consider  hydrophobic  effects  and  other  chemical 
interactions  between  molecules,  which  may  play  a  role  in 
defining  the  strength  properties.  This  is  a  limitation  of  the 
model;  however,  there  is  currently  no  method  to  explicitly 
include  these  effects  in  the  model  and  thus  this  task  is 
left  to  future  work.  We  expect,  however,  that  the  effect 
of  intermolecular  adhesion  is  limited  with  respect  to  the 
prediction  of  the  initial  strength  values.  This  is  based  on  our 
observation  that  in  AH-based  coiled-coil  protein  structures, 
failure  initiates  first  in  the  individual  AHs  (thus  defining  its 
strength  properties)  and  is  later  followed  by  uncoiling  of 
the  overall  coiled-coil  structure  [27].  This  suggests  that  the 
approach  taken  here  by  focusing  solely  on  H-bonds  is  a  good 
approximation  for  the  strength  properties  of  AH-based  protein 
domains. 


2.2.  Molecular  dynamics  simulation  procedure  and 
comparison  with  HBM 


4 


Atomistic  MD  with  explicit  solvent  is  used  to  simulate 
the  atomistic  mechanisms  that  control  the  deformation  and 
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Table  1.  Overview  of  important  variables  and  symbols  used  in  the 
manuscript. 


Mathematical 

symbol 

Unit 

Description 

kcal  moU^ 

Height  of  the  energy  barrier  at  the 
transition  state 

kcal  moU^ 

Energy  barrier  of  a  single  H-bond 

Vb 

A 

Location  of  the  energy  barrier 

/ 

pN 

Applied  force  at  molecule 

0 

°  (degrees) 

Angle  between  pulling  direction 
and  reaction  coordinate  of 
breaking  H-bond 

Icb 

JK-' 

Boltzmann  constant 

T 

K 

Absolute  temperature 

(jOq 

s-i 

Natural  bond  vibration  frequency 
(1  X  IQi^  s-i) 

X 

s“^ 

Off  rate:  bond  dissociation  per 
second 

V 

m  s“^ 

Macroscopically:  pulling  speed, 
microscopically  bond  breaking 
speed 

Vo 

m  s“^ 

Natural  bond  breaking  speed, 
when  no  load  is  applied 

£ 

% 

Engineering/molecular  strain: 
displacement  of  the  pulled  atom 
normalized  by  the  length  of  the 
molecule 

bi 

— 

Number  of  parallel  elements  at 
hierarchy  level  i 

ki 

Number  of  elements  at  hierarchy 
level  i  that  rupture 
simultaneously 

fv 

pN 

Force  contribution  as  a 
consequence  of  pulling  speed 

fhi 

pN 

Force  contribution  as  a 
consequence  of  hierarchy  level  i 

fhO 

pN 

Force  contribution  as  a 
consequence  of  the  basic 
hierarchy  (individual  H-bond) 

r 

% 

Robustness,  values  ranging  from 
0%  to  100% 

rupture  of  bonds  at  nanoscale,  and  are  used  here  to  validate 
the  predictions  put  forward  by  the  Hierarchical  Bell  model 
developed  above. 

MD  calculations  are  carried  out  using  NAMD  [28]  with 
a  CHARMM  [29]  force  field  with  explicit  water.  Each 
molecular  assembly  is  embedded  in  a  skin  of  explicit  TIPS 
water  at  pH  7.  The  single  AH  and  the  two-stranded  coiled- 
coils  (CC2)  are  taken  from  the  2B  segment  of  vimentin  IF 
(PDB  ID  1GK6,  residues  355-406).  The  four-stranded  CC 
(CC4)  is  taken  from  the  early  endosomal  SNARE  complex 
(PDB  ID  2NPS,  residues  A:  60-112,  B:  196-248,  C:  138- 
190,  D:  180-232).  Energy  minimization  is  carried  out  for 
100000  steps,  which  allows  for  a  favorable  conformation  to 
be  achieved.  This  is  followed  by  an  equilibration  procedure, 
during  which  the  molecular  assembly  is  heated  up  to  300  K 
with  a  rate  of  25  K  every  25  steps.  An  NVT  ensemble 
is  employed  to  hold  the  temperature  constant  at  the  final 
temperature  of  300  K.  Using  a  time  step  of  1  fs,  each  molecular 
is  subjected  to  equilibration  for  1  ns.  To  simulate  forced 
rupture  of  H-bonds  between  the  strands,  we  use  the  steered 
molecular  dynamics  procedure  (SMD)  with  a  constant  velocity 


pulling  scheme.  A  spring  constant  of  10  kcal  mol“^  A  ^  and 
a  pulling  rate  ranging  form  1  to  100  m  s“^  are  used  for  this 
purpose.  In  interatomic  potentials  (such  as  the  CHARMM 
force  field  as  used  for  the  studies  here),  H-bonds  are  modeled 
based  on  empirical  energy  expressions.  Such  expressions 
are  fitted  to  first  principles  or  experimental  measurements  of 
H-bond  geometries  and  energies.  In  CHARMM,  Coulomb 
and  Lennard- Jones  terms  between  the  donor  and  acceptor 
accurately  represent  the  H-bond  interaction,  making  an  explicit 
H-bond  interaction  formulation  unnecessary  [29].  Thereby, 
H-bond  models  in  such  force  fields  describe  the  ability  of 
donor-acceptor  pairs  of  H-bond  formation  to  interact  and 
form  intermolecular  or  intramolecular  bonding.  H-bonds  in 
CHARMM  force  fields  are  ‘reactive’,  that  is,  they  can  break 
and  reform  with  any  donor-acceptor  pair  in  a  simulation 
model.  This  resembles  the  conventional  way  of  thinking  of  H- 
bonds  as  weak  interactions  formed  within  and  between  protein 
domains. 

The  analysis  with  the  Hierarchical  Bell  model  shown  in 
figure  2(c)  is  carried  out  using  k\  =  \  according  to  the 
elementary  rupture  mode  under  the  loading  conditions  (see 
reference  [17]  for  details).  The  results  shown  in  this  article 
are  visualized  using  visual  molecular  dynamics  (VMD)  [30]. 

2.3.  Strength  analysis 

In  all  examples,  we  plot  the  force  per  AH,  that  is,  we 
normalize  the  strength  to  enable  a  suitable  comparison  between 
different  hierarchical  structures.  We  focus  on  the  hierarchical 
contribution  to  the  strength  and  robustness  of  the  protein 
structure  (the  contribution  to  strength  and  robustness  from 
the  pulling  velocity  has  the  same  absolute  value  at  a  fixed 
deformation  velocity,  for  different  structures). 

2.4.  Definition  of  robustness  in  biological  protein  constituents 

We  calculate  robustness  as  the  ratio  of  strength  of  a  defected 
system  and  an  intact  system,  by  following  Kitano’s  definition 
of  robustness  as  fault  insensitivity  [31].  For  the  calculations 
reported  below,  the  intact  system  is  defined  as  a  system  in 
which  all  AHs  contribute  to  strength,  whereas  in  the  defected 
system  all  except  one  AHs  (fault  on  second  hierarchy,  variation 
of  b2  [=number  of  AHs  in  a  bundle],  resulting  in  bf)  contribute 
to  the  strength.  The  robustness  is  defined  as  the  ratio  of  the 
strength  of  the  fiawed  system  divided  by  the  strength  of  the 
intact  system  (see  equation  (2)  for  the  definition  of  fhi), 

_  f{bi,ki,bl  =  b2-\) 

fho  +  fhiih  =3,br=3)  +  M{k2  =  \,b*  =  b2-l) 

N 

+  =  i,^i) 

i—3 
N 

+  =  hbi) 

/=2 

The  robustness  converges  towards  complete  fault  toler¬ 
ance  when  bi  oo,  as  shown  in  figure  4.  We  note  that  other 


II 
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Figure  3.  Performance  of  the  analyzed  structures  in  the  strength  per 
AH-robustness  space,  with  and  without  a  defect.  The  coarsely 
dashed  lines  represent  levels  of  equal  strength  (5)-robustness  (r) 
potential  (that  is,  the  product  of  both  values  is  equal  on  these  lines, 
r  •  s  =  const).  Robustness  and  strength  compete  on  these  lines.  The 
first  data  point  for  each  structure  represents  the  intact  system, 
whereas  the  second  data  point  shows  the  system  after  failure.  CC2s 
have  a  robustness  degree  of  80%  (robustness  equals  to  the  force  from 
hierarchical  strengthening  of  a  defect  system,  two  instead  of  three 
H-bonds  rupture  simultaneously,  divided  by  the  force  of  an  intact 
system,  when  all  three  H-bonds  rupture  at  once).  The  presence  of  a 
defect  moves  the  system  to  another  potential  line.  For  example,  due 
to  the  high  level  of  robustness,  the  CC4  structure  hardly  changes  its 
strength,  whereas  the  strength  of  a  single  AH  is  significantly  reduced. 
This  illustrates  how  robustness  and  strength  are  combined  in 
dependence  of  functional  requirements. 


robustness  perspectives  in  biological  materials  include  adap¬ 
tation  to  environmental  changes,  as  well  as  graceful  degrada¬ 
tion  [31].  This  definition  of  robustness  is  related  to  a  struc¬ 
ture’s  ability  to  tolerate  structural  defects;  it  is  therefore  related 
to  toughness.  Tough  materials  show  high  robustness  against 
catastrophic  failure  despite  the  presence  of  a  defect. 


3.  Computational  and  theoretical  results 

3.1.  Validation  of  the  Hierarchical  Bell  model  via  direct  MD 
simulation 

We  first  validate  the  Hierarchical  Bell  model  by  carrying  out 
full  atomistic  simulation  studies  in  explicit  water  of  the  rupture 
force  of  three  different  hierarchical  AH  protein  motifs  (see 
figure  2(a)  for  the  geometries).  Force-extension  histories  from 
the  MD  simulations  are  shown  in  figure  2(b).  Figure  2(c) 
shows  theoretical  predictions  of  the  strength  of  the  structures 
according  to  equation  (2)  (see  continuous  lines),  and  the 
direct  comparison  with  results  of  MD  simulations.  The 
predictions  based  on  the  model  require  only  two  fundamental 
input  parameters,  the  H-bond  energy  and  the  bond  breaking 
distance  Xb.  Notably,  these  values  are  the  same  for  all  structures 
shown  in  figure  2.  Thus,  the  analysis  shows  that  the  different 
strength  values  for  the  different  geometries  are  a  result  solely  of 
hierarchical  effects  resulting  from  different  geometries.  These 
simulations  confirm  that  the  formulation  of  the  Bell  model  for 
hierarchical  systems  given  in  equations  (2)  and  (3)  provides 
an  accurate  prediction  of  the  rupture  forces  or  strength  values 
based  on  the  particular  geometry  of  the  protein  structure. 


Strength  per  AH 


2  hiearchies 


Figure  4.  Robustness-strength  domain  for  different  hierarchical 
arrangements,  considering  the  effect  of  increasing  the  number  of 
subelements  on  the  particular  hierarchy  (arrow  points  in  direction  of 
increasing  number  of  hierarchies).  Subplot  (b)  depicts  examples  for 
2-hierarchy  and  3-hierarchy  systems,  each  with  one,  two  and  three 
elements.  These  results  show  that  in  general,  the  more  elements,  the 
more  robust  but  the  less  strong  (per  AH)  is  a  system.  This  also  shows 
the  tradeoff  between  robustness  and  strength.  The  strength  per  AH  is 
plotted,  to  enable  a  better  comparison  of  different  structures,  which 
feature  different  amounts  of  material  (that  is,  number  of  AHs)  per 
cross-sectional  area.  This  equals  to  a  normalization  of  force  by  the 
cross-sectional  area,  which  leads  to  the  strength  of  a  material. 


3.2.  Applications  and  case  studies 

The  Hierarchical  Bell  model  is  now  used  in  a  series  of 
theoretical  analyses  of  three  case  studies,  addressing  the 
question,  how  are  AH-based  protein  domains  capable  of 
unifying  strength  and  robustness?  The  direct  atomistic 
simulation  of  the  systems  considered  here  is  not  possible  due 
to  computational  limitations,  so  that  the  analysis  is  limited  to 
theoretical  considerations. 

For  the  following  analysis  we  use  the  from  the  single 
AH  (E"^  =  4.21  kcal  mol“^  and  Vb  =  1.21  A)  as  extracted 
from  full  atomistic  simulations  of  AH  structures  [17].  These 
values  are  close  to  results  of  experimental  and  theoretical 
investigations  of  the  rupture  mechanism  of  AH  protein 
domains  [17,  32].  Earlier  full  atomistic  studies  in  explicit 
water  solvent  by  the  same  authors  have  shown  that  the 
fundamental  fracture  mechanism  of  a  single  AH  element  is 
the  simultaneous  rupture  of  approximately  3  H-bonds  [17]. 
This  mode  represents  the  basic  unit  mechanism  of  rupture  in 
vivo,  at  pulling  rates  below  0.1  m  s“^  [17].  This  mechanism 
is  considered  in  the  case  studies  reported  here  (figures  3-5), 
represented  by  =  3,  thereby  providing  a  realistic  description 
of  the  behavior  closer  to  experimentally  accessible  pulling 
rates  (it  was  shown  in  an  earlier  analysis  [17]  that  this  mode 
of  deformation,  k\  =3,  dominates  at  low  deformation  rates, 
instead  of  =  1,  which  is  due  to  the  high  pulling  rates  in  MD 
simulations). 

We  now  apply  the  Hierarchical  Bell  model  to  three  AH- 
based  structures  (these  are  the  same  structures  as  discussed  in 
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Figure  5.  This  figure  gives  an  example  of  four  different  structures 
with  the  same  number  of  subelements  (that  is,  eight  AHs)  but  in 
different  hierarchical  arrangements.  Subplot  (a)  shows  the  four 
different  architectures.  For  simplicity,  individual  H-bonds  on  the 
lowest  hierarchical  scale  are  not  shown;  instead  one  line  represents 
three  H-bonds  as  one  AH.  Subplot  (b)  shows  the  concurrence 
between  strength  and  robustness,  which  depends  on  the  degree  of 
redundancies  on  the  different  hierarchical  levels.  The  level  of 
robustness  increases  with  increasing  redundancies  on  a  particular 
level.  Dependent  on  the  hierarchical  arrangement  of  the  elements, 
different  potentials  of  strength  and  robustness  can  be  reached. 
Subplot  (c)  shows  the  contributions  of  each  hierarchy  to  the  overall 
strength  (not  strength  per  AH,  as  shown  in  subplot  (b)).  As  we 
assume  that  in  each  AH  three  H-bonds  break  simultaneously,  each 
structure  (featuring  an  AH  as  the  smallest  subelement)  has  the  same 
contribution  from  hierarchy  0.  This  is  also  the  highest  amount  of 
strength  contribution  and  shows  the  significance  of  the  strength  of 
H-bonds,  which  depends  on  the  solvent  and  the  environment.  The 
other  contributions  are  of  hierarchical  origin.  The  force  contribution 
from  hierarchy  1  is  zero,  since  3  out  of  3  H-bonds  break,  which 
lowers  the  logarithmic  multiplicator  to  zero. 


the  validation  of  the  model,  shown  in  figure  2(a)).  Figure  3 
shows  the  performance  of  the  analyzed  structures  in  the 
strength  per  AH-robustness  space  with  and  without  a  defect 
of  one  basic  element  (see  section  2.4  for  the  definition  of 
robustness;  typical  strength  values  under  slow  deformation 
conditions  [17]  range  between  100  and  several  hundred  pN, 
at  deformation  speeds  below  0.1  m  s“^).  A  system  with  a  high 
level  of  robustness  and  thus  with  a  high  level  of  redundancies 


hardly  changes  the  strength  even  when  a  defect  appears  (as 
seen  in  the  4- stranded  CC,  CC4).  In  contrast,  a  system  with 
a  low  level  of  robustness  and  thus  less  redundancies  (as  seen 
in  the  single  AH)  has  a  significantly  reduced  strength  when 
a  defect  appears.  However,  increasing  robustness  goes  along 
with  a  loss  of  strength  per  element  (that  is,  per  AH),  since  the 
overall  strength  of  a  structure  is  not  directly  proportional  to 
the  number  of  parallel  elements.  Theoretically,  a  robustness 
value  of  100% — leading  to  the  highest  strength  of  a  defected 
system — would  appear  for  an  infinite  number  of  elements  at 
each  hierarchical  scale. 

Figure  4  depicts  a  systematic  investigation  of  the 
robustness-strength  behavior  under  an  increasing  number  of 
elements  at  each  hierarchy,  illustrating  this  effect.  Each  line 
represents  one  level  of  hierarchies  (e.g.  2  hierarchies  equal 
to  the  AH  level,  3  hierarchies  equal  to  the  CC  level),  where 
the  number  of  elements  on  this  particular  hierarchies  is  varied 
(e.g.  one  AH,  two  AHs,  for  2  hierarchies  system  and  one  CC, 
two  CCs,  and  others,  for  a  three-hierarchy  system).  We  find 
that  at  each  equihierarchical  line,  with  an  increasing  number  of 
elements,  the  strength  per  AH  decreases,  while  the  robustness 
increases.  However,  increasing  the  number  of  elements  at  a 
specific  hierarchical  scale  is  inefficient,  as  it  leads  to  extensive 
material  use  and  a  decrease  in  the  overall  strength.  Therefore, 
the  introduction  of  hierarchies,  enabling  the  optimization  of 
strength  and  robustness  under  limited  material  use  becomes 
significant. 

To  illustrate  this  point,  we  arrange  eight  single  AHs  in 
different  hierarchical  structures — asking  the  question:  how  can 
one  arrange  eight  AHs  to  obtain  different  levels  of  robustness 
and  strength?  As  shown  in  figure  5,  the  systems  (schematics 
in  figure  5(a))  consist  of  two,  three  and  four  hierarchies. 
The  differences  in  robustness  and  strength  calculated  with  the 
Hierarchical  Bell  model  are  not  achieved  through  additional 
use  of  materials,  but  purely  through  different  hierarchical 
arrangements  (figures  5(b)  and  (c)).  To  the  best  of  our 
knowledge,  this  tuning  of  properties  in  the  strength-robustness 
domain  as  illustrated  here  has  been  shown  theoretically  here 
for  the  first  time.  In  the  robustness-strength  map,  the  ‘best’ 
material  behavior  is  the  one  in  which  high  robustness  is 
achieved  at  large  strength — referred  to  as  a  ‘high  potential’.  It 
can  be  seen  that  system  2  has  the  highest  potential.  Notably,  it 
is  not  the  system  with  the  highest  hierarchical  level  (system  4), 
nor  the  system  with  the  highest  level  of  redundancies  (system 
1).  In  other  words,  system  2  features  the  best  combination 
of  redundancies  at  different  hierarchical  levels,  and  features 
superior  mechanical  performance  without  additional  material 
use. 

The  analysis  reported  here  shows  that  with  different 
arrangements  (by  changing  the  number  of  subelements  as 
shown  in  figure  4,  and  by  changing  the  hierarchical  geometries 
as  shown  in  figure  5),  almost  any  point  in  the  strength- 
robustness  space  can  be  achieved.  These  results  illustrate 
how  AH  protein  domains  solve  the  conflict  between  strength 
and  robustness,  by  introducing  hierarchies  as  an  additional 
design  variable.  This  finding  is  the  most  important  result 
of  the  study  reported  in  this  paper.  The  results  shown  in 
figure  5  results  suggest  that  the  level  of  hierarchical  depth 
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and  strength  may  be  balanced  in  biological  protein  materials, 
allowing  biological  materials  to  maximize  the  mechanical 
performance  while  minimizing  the  use  of  materials.  In 
agreement  with  this  notion,  proteomics  analysis  reveals  that 
CC2  structures  (see  figure  2(a))  are  most  common  in  biology 
(due  their  dominance  in  intermediate  filaments,  feather,  hair, 
hoof  and  other  materials  [7,  8, 11]),  maybe  since  they  provide  a 
compromise  between  strength  and  robustness  (having  two  AHs 
instead  of  one  is  the  minimum  level  of  possible  redundancies). 
Further,  the  energy  dissipation  per  AH  in  the  plateau  regime  II 
(see,  e.g.  figure  2(b))  depends  directly  on  the  force  level  and 
might  be  another  reason  for  the  dominance  of  CC2  structures, 
as  they  feature  a  minimum  level  of  redundancies  at  high  levels 
of  energy  dissipation  density. 

4.  Discussion  and  conclusions 

We  have  shown  here  that  with  different  structural  arrange¬ 
ments,  different  combinations  of  strength  and  robustness  can 
be  achieved.  This  finding  is  the  most  important  result  of  the 
case  studies  put  forth  in  this  article:  it  illustrates  that  the  con¬ 
flict  between  strength  and  robustness  can  be  resolved  by  intro¬ 
ducing  hierarchies  as  an  additional  design  variable.  This  pro¬ 
vides  important  insight  into  structure-property  relationships  in 
protein  materials,  contributing  to  on-going  efforts  at  the  inter¬ 
face  of  materials  science  and  biology  [33].  These  results  fur¬ 
ther  suggest  that  the  level  of  hierarchical  depth  and  strength 
may  be  balanced  in  biological  protein  materials,  since  the 
robustness  and  strength  are  not  completely  inversely  propor¬ 
tional,  allowing  biological  materials  to  maximize  the  mechan¬ 
ical  performance  while  minimizing  the  use  of  materials.  Over¬ 
all  our  analysis  illustrates  that  the  introduction  of  hierarchies 
is  the  key  to  unify  disparate  material  properties.  Applying 
this  insight  to  the  design  of  materials  will  allow  an  extended 
use  of  hierarchies  in  bioinspired  or  biomimetic  synthetic  ma¬ 
terials  at  nanoscale,  such  as  hierarchically  organized  CNT- 
bundles,  nano  wires,  CNT-protein  or  polymer-protein  compos¬ 
ites  [34-36].  The  combination  of  synthetic  and  natural  consti¬ 
tuting  elements  (e.g.  proteins)  could  be  a  particularly  promis¬ 
ing  strategy. 

The  increasing  use  of  protein  building  blocks  in  the 
development  of  novel  nanomaterials,  such  as  nano  wires, 
nanotubes,  and  others  [13]  requires  the  development  of  new 
engineering  models  that  enable  the  systematic  design  of  the 
use  of  nanoscale  constituents  in  the  makeup  of  larger-scale 
materials.  Our  model  provides  such  insight  and  design  rules, 
here  shown  for  the  example  of  AH  structures,  which  facilitate 
the  development  of  novel  nanostructures  based  on  proteins. 
Most  importantly,  our  model  explains  a  fundamentally  new 
concept,  that  is,  by  simply  rearranging  the  same  number 
of  nanoscale  elements  into  hierarchies,  one  can  change  the 
performance  of  the  material  in  the  strength-robustness  space 
(see  figure  5).  This  makes  the  continuous  invention  of  new 
basic  building  blocks  unnecessary.  Thus  the  broad  application 
of  universal  building  blocks  in  highly  diverse  architectures 
might  be  a  biological  strategy  that  enables  adaptation  to 
changes  in  the  environment  directly  by  adopting  the  structural 
arrangement  of  the  same  basic  building  blocks.  This  concept 


also  appears  at  the  level  of  primary  structure  of  proteins. 
Instead  of  inventing  new  amino  acids,  a  limited  number  of 
universal  20  amino  acids  are  combined  in  different  ways 
to  create  functional  complexity,  accounting  for  the  great 
variety  of  biological  protein  materials.  These  types  of 
concepts  represent  an  opportunity  for  further  studies  that  could 
investigate  this  hypothesis  in  greater  detail.  In  additional  to 
the  biological  context,  detailed  studies  of  this  concept  as  a 
novel  ‘engineering  paradigm’  are  crucial  to  advance  this  field. 
This  might  result  in  synthetic  materials  as  mechanomutable, 
smart  structures,  which  continuously  and  independently  adapt 
to  environmental  changes  at  each  length-  and  timescale, 
consisting  of  a  handful  of  ‘universal’  building  blocks  at  the 
nanoscale. 

Even  though  our  model  is  focused  on  AH  protein 
structures,  the  main  results  of  the  work  should  be  generally 
valid  for  other  protein  structures  that  are  primarily  stabilized 
by  H-bonds,  in  particular  the  finding  that  hierarchical 
arrangements  of  H-bonds  are  crucial  for  the  ability  of  protein 
materials  to  combine  strength  and  robustness.  This  property 
of  biological  protein  materials  has  often  been  pointed  out  in 
the  literature  [2,  37],  but  has  not  yet  been  explained  based  on 
a  fundamental  physics  model  as  accomplished  here  with  the 
Hierarchical  Bell  model  (equations  (2)  and  (3)).  The  observed 
behavior  provides  further  evidence  that  this  behavior  may  have 
its  origin  at  the  molecular  scale.  This  insight  could  be  crucial  to 
translate  biological  material  concepts  towards  nanotechnology 
applications.  However,  caution  must  be  taken  since  other 
types  of  chemical  bonding  as  well  as  the  environmental 
conditions  such  as  temperature  or  pH  (possibly  strengthening 
and  weakening  different  types  of  bonds)  may  also  contribute 
to  a  protein’s  mechanical  stability  (see  also  discussion  in  the 
introduction),  and  that  these  models  should  be  considered 
appropriately.  Further  studies,  perhaps  using  full  atomistic 
MD  with  different  model  protein  structures  in  explicit  solvent, 
could  provide  crucial  insight  into  these  aspects. 

Figure  6  shows  how  our  theory  can  be  applied  in  an 
analysis  of  the  mechanical  performance  of  a  different  protein 
nanostructure,  a  triple  beta-helix  amyloid  protein  nanotube 
that  is  part  of  the  T4  bacteriophage  virus  (structure  shown  in 
figure  6(a)).  This  amyloid  structure  is  able  to  withstanding 
large  pulling  and  compression  forces  in  the  order  of  several  nN 
[38],  and  might  be  of  great  interest  for  future  mechanical  and 
nanoelectronic  applications  [13,  39].  Each  convolution  in  this 
structure  consists  of  three  beta- sheets  arranged  in  an  equilateral 
triangle,  where  each  beta- sheet  consists  of  a  cluster  of  several 
H-bonds.  The  structure  therefore  consists  of  two  hierarchies 
per  convolution.  Figure  6(b)  depicts  the  representation  of  this 
protein  structure  in  the  Hierarchical  Bell  model.  How  would 
the  performance  of  this  structure  change  if  the  structure  is 
altered?  Figure  6(c)  depicts  a  representation  of  a  system  that 
would  contain  all  H-bonds  per  convolution  in  one  large  cluster, 
that  is,  this  structure  has  lost  the  arrangement  of  H-bonds  in 
clusters.  Applying  equation  (3),  and  =  5  kcal  mol“\ 
and  assuming  that  the  bonds  rupture  one  by  one  (ki  =  1 
as  expected  from  MD  simulations  of  similar  structures),  the 
Hierarchical  Bell  model  predicts  that  the  introduction  of  the 
hierarchical  level  by  clustering  H-bonds  increases  the  effective 
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1  beta-sheet  1  Convolution  =  3  Sheets 


1  Convolution 

Figure  6.  The  applicability  of  our  model  to  a  different  protein  nanostructure,  a  triple  beta-helix,  as  it  could  be  used  for  mechanical  and 
electronic  nanodevices.  The  analysis  shown  here  is  for  the  needle  of  a  T4  bacteriophage.  Subplot  (a)  depicts  the  protein’s  ribbon  structure 
from  two  different  views.  Each  convolution  consists  of  three  beta- sheets  arranged  in  an  equilateral  triangle,  where  each  beta- sheet  consists  of 
a  cluster  of  several  H-bonds.  Subplot  (b)  depicts  the  representation  of  this  protein  structure.  The  structure  consists  of  two  hierarchies  per 
convolution,  where  H-bonds  are  arranged  in  three  clusters  in  each  convolution.  Subplot  (c)  depicts  the  representation  of  a  system  that  would 
contain  nine  parallel  bonds  per  convolution,  that  is,  all  available  H-bonds  per  convolution  are  arranged  in  one  large  cluster.  Our  model  predicts 
that  the  height  of  the  effective  energy  barrier  Eb  of  system  (b)  is  28%  higher  than  the  one  of  system  (c),  illustrating  the  effect  of  using 
hierarchies  to  improve  the  performance  of  nanostructures. 


energy  barrier  E'b  by  28%  compared  to  the  arrangement 
with  only  one  hierarchy,  leading  to  greater  chemomechanical 
stability.  This  illustrates  the  effect  of  using  hierarchies  to 
improve  the  performance  of  nanostructures.  This  simple 
analysis  illustrates  for  a  specific  example  how  our  model  could 
be  used  in  the  design  of  novel  protein-based  nanomaterials. 

The  failure  of  engineering  materials  and  structures  has 
been  studied  extensively  and  has  impacted  our  world  by 
enabling  the  design  of  complex  structures  such  as  buildings, 
airplanes,  cars  and  devices.  However,  the  mechanisms 
of  failure  in  biological  materials  and  how  it  leads  to 
the  breakdown  of  components  in  our  body  is  not  well 
understood.  Thus,  characterizing  how  protein  materials 
fail  has  significant  implications  that  may  eventually  lead  to 
an  improved  understanding  of  diseases  and  injuries.  The 
link  between  molecular  structure  and  material  properties  in 
the  strength-robustness  domain  may  lead  to  a  paradigm 
shift  in  the  understanding  of  which  physical  mechanisms 
govern  the  behavior  of  biological  systems.  Understanding 
the  fundamental  physical  laws  that  control  the  properties 
of  hierarchical  protein  materials  enables  us  also  to  link  the 
structural  protein  organization  to  the  appropriate  biological 
functions.  Engineering  materials  using  a  bottom-up  approach 
that  begins  at  the  atomistic  level,  inspired  by  biological  protein 
material  concepts,  may  transcend  the  borders  that  currently  lie 
between  life  sciences  and  engineering.  The  transfer  towards 
the  design  of  novel  nanostructures  facilitates  the  development 
of  de  novo  multi-functional  and  mechanically  active,  tunable 


and  changeable  materials  [1,  40-44],  for  example  new 
organic  and  organic-inorganic  composites  that  primarily 
consist  of  chemical  elements  that  appear  in  our  environment 
in  practically  unlimited  amount  (chemical  elements  such  as  C, 
H,  N,  O,  S). 

This  may  lay  the  foundation  for  a  new  engineering 
paradigm  that  includes  the  design  of  structures  and  materials 
starting  at  the  molecular  level,  from  bottom-up,  to  the 
macroscale,  to  create  new  materials  and  structures  that  mimic 
and  exceed  the  properties  found  in  biological  analogs.  The 
development  of  a  fundamental  science  driven  framework  that 
involves  a  solid  understanding  of  fundamental  concepts  is 
crucial  for  studies  of  biological  systems,  disease  diagnosis 
and  treatment,  as  well  as  the  design  of  novel  biomaterials. 
It  is  the  key  to  reverse-engineer  the  human  body,  the  key 
to  understand  diseases  at  multi- scale  levels  (cancer,  genetic 
diseases,  infectious  diseases),  to  enable  advanced  treatments 
(intervention  at  level  of  relevant  proteins,  nanomedicine, 
protein  hierarchies,  biomaterials)  and  diagnostics  (mechanical 
disease  signature,  protein  misregulation). 

Recently  it  was  reported  that  the  first  bundles  of 
coiled-coil  proteins  were  generated  synthetically  [40],  and 
many  other  peptide  synthesis  techniques  are  progressing 
rapidly  [1,  40,  42-44],  allowing  to  create  self-assembled 
nanostructures  of  amyloid  fibers  [13]  or  collagen  fibrils  [45]. 
The  model  reported  in  this  paper,  combined  with  these  new 
manufacturing  techniques,  may  be  the  first  step  towards  a  de 
novo  bottom-up  structural  design  of  protein-based  materials.  In 
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addition  to  protein  materials  this  theory  could  be  applicable  to 
other  nanoscale  devices,  which  exhibit  a  hierarchical  structure 
and  are  governed  by  stochastic  processes  of  failure,  such  as 
polymer  brushes  or  multi-layer  films.  A  detailed  analysis  of  the 
applicability  of  this  theory  to  these  structures  is  left  to  future 
work. 

The  field  of  genomics  is  concerned  with  the  study  of 
genes  and  their  effects  on  macroscopic  functions,  and  has  led 
to  considerable  medical  advances.  Genomics,  however,  does 
not  elucidate  material  properties,  nor  the  mechanistic  relation 
of  hierarchical  multi- scale  structures  and  their  resulting 
properties.  The  multi-scale  behavior  of  protein  assemblies 
with  the  goal  of  elucidating  the  relation  between  structure 
and  material  properties  represents  a  grand  challenge  at  the 
interface  of  materials  science  and  biology.  This  gap  in 
understanding  could  be  closed  by  systematically  studying  the 
material  properties  of  hierarchical  protein  structures,  their 
effect  on  the  macroscopic  properties  (through  development 
of  structure-property  relationships),  and  the  role  of  material 
properties  in  their  biological  context,  an  effort  defined  as 
materiomics  [19]. 

Acknowledgments 

This  research  was  supported  by  the  Army  Research  Office, 
grant  #  W911NF-06- 1-0291  (program  officer  Dr  Bruce 
LaMattina),  by  a  National  Science  Foundation  CAREER 
Award  (program  manager  Dr  Jimmy  Hsia),  grant  #  0642545, 
and  by  a  grant  from  the  Air  Force  Office  of  Scientific 
Research,  grant  #  FA9550-08- 1-0321  (program  manager 
Dr  Les  Lee).  TA  acknowledges  support  from  the 
German  National  Academic  Foundation  (Studienstiftung  des 
Deutschen  Volkes),  the  Hamburg  Foundation  for  research 
studies  abroad  (Hamburger  Stipendienprogramm)  and  the  Dr. 
Juergen  Ulderup  Foundation.  The  authors  thank  Professor 
Lothar  Gaul  (University  Stuttgart)  and  Professor  Reinhard 
Lipowsky  (Max-Planck  Institute  of  Colloids  and  Interfaces 
Potsdam)  for  their  continuous  interest  and  support  of  our  work. 

References 

[1]  Langer  R  and  Tirrell  D  A  2004  Designing  materials  for  biology 

and  medicine  Nature  428  487-92 

[2]  Fratzl  P  and  Weinkamer  R  2007  Nature’s  hierarchical  materials 

Prog.  Mater.  Sci.  52  1263-334 

[3]  Ingber  D  E  2006  Cellular  mechanotransduction:  putting  all  the 

pieces  together  again  FASEB  J.  20  811-27 

[4]  Hirth  J  P  and  Lothe  J  1982  Theory  of  Dislocations  (New  York: 

Wiley-Interscience) 

[5]  Murphy  K  P  2001  Protein  Structure,  Stability,  and  Folding 

{Methods  in  Molecular  Biology  vol  168)  ed  N  W  John 

(Totowa:  Humana  Press) 

[6]  Gruber  M  and  Lupas  A  N  2003  Historical  review:  Another  50th 

anniversary-new  periodicities  in  coiled  coils  Trends 

Biochem.  Sci.  28  679-85 

[7]  Herrmann  H  and  Aebi  U  2004  Intermediate  filaments: 

molecular  structure,  assembly  mechanism,  and  integration 

into  functionally  distinct  intracellular  scaffolds  Annu.  Rev. 

Biochem.  73  749-89 

[8]  Alberts  B  et  al  2002  Molecular  Biology  of  the  Cell  (London: 

Taylor  and  Francis) 


[9]  Keten  S  and  Buehler  M  J  2008  Asymptotic  strength  limit  of 
hydrogen-bond  assemblies  in  proteins  at  vanishing  pulling 
rates  Phys.  Rev.  Lett.  100  198301 

[10]  Lang  M  2008  Lighting  up  the  mechanome  Frontiers  of 

Engineering:  Reports  on  Leading-Edge  Engineering  from 
the  2007  Symp.  National  Academy  of  Engineering  of  the 
National  Academies 

[11]  Rowat  AC  et  al  2008  Towards  an  integrated  understanding  of 

the  structure  and  mechanics  of  the  cell  nucleus  Bioessays 
30  226-36 

[12]  Ackbarow  T  and  Buehler  M  J  2008  Hierarchical  coexistence  of 

universality  and  diversity  controls  robustness  and 
multi-functionality  in  protein  materials  Theor.  Comput. 

Nano  sci.  5  1193-204 

[13]  Scheibel  T  et  al  2003  Conducting  nanowires  built  by  controlled 

self-assembly  of  amyloid  fibers  and  selective  metal 
deposition  Proc.  Natl  Acad.  Sci.  USA  100  4527-32 

[14]  Ahl  V  and  Allen  TEH  1996  Hierarchy  Theory -A  Vision, 

Vocabulary,  and  Epistemology  (New  York:  Columbia 
University  Press) 

[15]  Buehler  M  J  and  Ackbarow  T  2007  Fracture  mechanics  of 

protein  materials  Mater.  Today  10  46-58 

[16]  Buehler  M  J  and  Gao  H  2006  Dynamical  fracture  instabilities 

due  to  local  hyperelasticity  at  crack  tips  Nature  439  307-10 

[17]  Ackbarow  T  et  al  2007  Hierarchies,  multiple  energy  barriers 

and  robustness  govern  the  fracture  mechanics  of 
alpha-helical  and  beta-sheet  protein  domains  Proc.  Natl. 
Acad.  Sci.  USA  104  16410-5 

[18]  Keten  S  and  Buehler  M  J  2008  Geometric  confinement  governs 

the  rupture  strength  of  H-bond  assemblies  at  a  critical  length 
scale  Nano  Lett.  8  743-8 

[19]  Buehler  M  J  and  Keten  S  2008  Elasticity,  strength  and 

resilience:  A  comparative  study  on  mechanical  signatures  of 
O' -helix,  yS- sheet  and  tropocollagen  domains  Nano  Res. 

1  63-71 

[20]  Marszalek  VEet  al  1999  Mechanical  unfolding  intermediates 

in  titin  modules  Nature  402  100-3 

[21]  Evans  E  and  Ritchie  K  1997  Dynamic  strength  of  molecular 

adhesion  bonds  J.  72  1541-55 

[22]  Sheu  S-Y  et  al  2003  Energetics  of  hydrogen  bonds  in  peptides 

Proc.  Natl  Acad.  Sci.  100  12683-7 

[23]  Seifert  U  2000  Rupture  of  multiple  parallel  molecular  bonds 

under  dynamic  loading  Phys.  Rev.  Lett.  84  2750-3 

[24]  Dudko  O  K,  Hummer  G  and  Szabo  A  2006  Intrinsic  rates  and 

activation  free  energies  from  single-molecule  pulling 
experiments  P/zy 5.  Rev.  Lett.  96  108101 

[25]  Sotomayor  M  and  Schulten  K  2007  Single-molecule 

experiments  in  vitro  and  in  silico  Science  316  1 144-8 

[26]  Bell  G 1  1978  Models  for  specific  adhesion  of  cells  to  cells 

Science  200  618-27 

[27]  Ackbarow  T  and  Buehler  M  J  2007  Superelasticity,  energy 

dissipation  and  strain  hardening  of  vimentin  coiled-coil 
intermediate  filaments:  Atomistic  and  continuum  studies 
J.  Mater.  Sci.  42  8771-87 

[28]  Nelson  MT  et  al  1996  NAMD:  a  parallel,  object  oriented 

molecular  dynamics  program  Int.  J.  Supercomput.  Appl. 

High  Perform.  Comput.  10  251-68 

[29]  MacKerell  AD  et  al  1998  All-atom  empirical  potential  for 

molecular  modeling  and  dynamics  studies  of  proteins 
J.  Phys.  Chem.  B  102  3586-616 

[30]  Humphrey  W,  Dalke  A  and  Schulten  K  1996  VMD:  visual 

molecular  dynamics  J.  Mole.  Graph.  14  33 

[31]  Kitano  H  2002  Computational  systems  biology  Nature 

420  206-10 

[32]  Sheu  SY  et  al  2003  Energetics  of  hydrogen  bonds  in  peptides 

Proc.  Natl  Acad.  Sci.  USA  100  12683-7 

[33]  Ackbarow  T  and  Buehler  M  J  2008  Nanomechanical  strength 

mechanisms  of  hierarchical  biological  materials  and  tissues 
Comp.  Methods  Biomech.  Biomed.  Eng.  11  595-607 


10 


Nanotechnology  20(2009)075103 


T  Ackbarow  and  M  J  Buehler 


[34]  Cui  XQetal  2007  Biocatalytic  generation  of 

ppy-enzyme-CNT  nanocomposite:  from  network  assembly 
to  film  growth  J.  Phys.  Chem.  C  111  2025-31 

[35]  Hule  R  A  and  Pochan  D  J  2007  Polymer  nanocomposites  for 

biomedical  application  MRS  Bull  32  5 

[36]  Winey  K I  and  Vaia  R  A  2007  Polymer  nanocomposites  MRS 

Bull  32  5 

[37]  Fratzl  P  et  al  2004  Structure  and  mechanical  quality  of  the 

collagen-mineral  nano-composite  in  bone  J.  Mater.  Chem. 
14  2115-23 

[38]  Keten  S  and  Buehler  M  J  2008  Large  deformation  and  fracture 

mechanics  of  a  beta-helical  protein  nanotube:  atomistic  and 
continuum  modeling  Comput.  Methods  Appl.  Mech.  Eng. 
197  3203-14 


[39]  Barrau  S  et  al  2008  Integration  of  amyloid  nanowires  in  organic 

solar  cells  Appl.  Phys.  Lett.  93  023307 

[40]  Papapostolou  D  et  al  2007  Engineering  nanoscale  order  into  a 

designed  protein  fiber  Proc.  Natl  Acad.  Sci.  0700801 104 

[41]  Petka  W  A  et  al  1998  Reversible  hydrogels  from 

self-assembling  artificial  proteins  Science  281  389-92 

[42]  Smeenk  ]  Metal  2005  Controlled  assembly  of  macromolecular 

beta- sheet  fibrils  Ang^w.  Chem.  Int.  Edn  44  1968-71 

[43]  Zhao  X  J  and  Zhang  S  G  2006  Molecular  designer 

self-assembling  peptides  Chem.  Soc.  Rev.  35  1105-10 

[44]  Mershin  A  et  al  2005  A  classic  assembly  of  nanobiomaterials 

Nat.  Biotechnol.  23  1379-80 

[45]  Kotch  F  W  and  Raines  R  T  2006  Self-assembly  of  synthetic 

collagen  triple  helices  Proc.  Natl  Acad.  Sci.  USA 
103  3028-33 


11 


